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ABSTRACT 



This thesis considers the influence of four secondary effects 
on the analysis of a simple ell-shaped piping structure when it is 
uniformly heated. The four effects are: (a) axial deformation, 

(b) shearing deformation, (c) beam-column effect, and (d) dif- 
ference between arc and chord. They are usually neglected in a 
piping flexibility analysis. 

In a paper of several years ago, J. E. Brock developed a theory 



conclusions. The present thesis reviews Brock* s theory and corrects 
some errors in coefficients. It then describes a digital computer 
program for IBM System 360 which performs the corresponding calculations 
Finally, it draws the new and quite significant conclusion that con- 
ventional piping stress analysis, which neglects these secondary effects 
can result in gross errors in evaluating stresses in piping config- 



for this investigation but di 






t any numerical results or 



urations which are likely to be encountered in piping practice. 



NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIF. 93940 



TABLE OF CONTENTS 



Section Page 

I INTRODUCTION 9 

1-1 Background 9 

1-2 Objectives of Work 10 

1-3 Assumptions and Limitations 10 

1- 4 Notation 11 

II THEORETICAL ANALYSIS AND INVESTIGATION TECHNIQUE 14 

2- 1 Review of Flexibility Analysis for the Ell- 14 

Shaped Pipe 

2-2 Dimensionless Parameters 19 



2-3 Requirements Imposed on the Digital Computer 20 



Program 

III DESCRIPTION OF COMPUTER PROGRAM 21 

3-1 General Remarks 21 

3-2 Input 22 

3-3 Loading of Parameter Values 28 

3-4 Combination of Secondary Effects 29 

3-5 Computation by Iteration 30 

3-6 Incrementing the Parameter 32 

3-7 Output 34 

3- 8 Subroutines 38 

IV RESULTS OF INVESTIGATION 41 

4- 1 Errors in ASME Paper No. 63 - APMW - 21 41 

4-2 Numerical Calculations 42 

4-3 Final Results and Conclusions 42 



3 



TABLE OF CONTENTS (cont.) 



Section Page 

BIBLIOGRAPHY 52 

APPENDIX 

A COMPUTER PROGRAM NOMENCLATURE 

A- I Main Program Nomenclature 53 

A-2 Subroutine ALPHA Nomenclature 57 

A- 3 Subroutine FACTOR Nomenclature 57 

A-4 Subroutine DRAW Description 58 

B COMPUTER PROGRAM FLOW CHART 65 

C DATA CARD DESCRIPTION 67 

C-l Data Card Placement 67 

C-2 Sample Data Card Set 69 

D CONTROL CARD ARRANGEMENT 70 

E PROGRAM LISTING 71 

E-l Main Program 71 

E-2 Subroutine ALPHA 84 

E-3 Subroutine FACTOR 85 

E-4 Sample Output 88 

F SAMPLE CALCULATION 89 

G ARC AND CHORD DIFFERENCE 91 



4 



LIST OF TABLES 



Page 



Table I Notation 11 

Table II Formulas for the coefficients 13 

Table III Expression relating the quantities of 17 

general case with those of the ell-shaped 
pipe 

Table IV Relates Factor to its secondary effect 39 

Table V List of combination designators with 40 

associated values of MT & Factor Values 



5 



LIST OF ILLUSTRATIONS (FIGURES) 



Page 



Figure 


1 


Uniform Tip Loaded Cantilever 


14 


Figure 


2 


Ell- shaped pipe experiencing no 
deformation due to temperature 
difference 


16 


Figure 


3 


Ell-shaped pipe experiencing 
deformation due to temperature 
difference 


16 


Figure 


4 


Free body diagram of both legs of the 
ell- shaped pipe under heating 


17 


Figure 


5 


Illustration showing the two different 
parameter sizes 


34 


Figure 


6 


Moment/Exact Moment vs. leg ratio for 
Arc and Chord only suppressed 


44 


Figure 


7 


Moment/Exact Moment vs. leg ratio for 
Arc and Chord only allowed 


45 


Figure 


8 


Moment /Exact Moment vs. leg ratio for 
Shear only suppressed 


46 


Figure 


9 


Moment/Exact Moment vs. leg ratio for 
Shear only allowed 


47 


Figure 


10 


Moment /Exact Moment vs. leg ratio for 
Axial deformation only suppressed 


48 


Figure 


11 


Moment/Exact Moment vs. leg ratio for 
Axial deformation only allowed 


49 


Figure 


12 


Moment/Exact Moment vs. leg ratio for 
Beam column only suppressed 


50 


Figure 


13 


Moment/Exact Moment vs. leg ratio for 
Beam column only allowed 


51 


Figure 


G-l 


Arc and Chord Difference 


91 



7 



ACKNOWLEDGEMENTS 



A large degree of credit for this thesis is due to 
Professor John E. Brock of the Mechanical Engineering Department 
of the Naval Postgraduate School. Some of his work provided the 
source material for the work reported here and his direction, 
assistance, and encouragement during the course of this work were 
contributions invaluable to the author. 

The author expresses most sincere appreciation to 
Professor John E. Brock. 



8 



SECTION I 



INTRODUCTION 



1-1 Background 

In 1964 Professor John E. Brock of the Mechanical Engineering 
Department of the Naval Postgraduate School, Monterey, California, 
published a paper entitled, M Some Secondary Effects in a Simple 
Piping Structure Under Heating".^ 

The secondary effects referred to are axial deformation, shear, 
beam column action, and difference between arc and chord. The 
structural element considered is an ell-shaped piping configuration 
with a square corner. 

In the normal analysis of such a structure, the above mentioned 
secondary effects are neglected, being regarded as insignificant in 
relieving the stresses caused by thermal expansion. Only the 
principal sources of relief, bending and torsion, are usually consid- 
ered. 

Professor Brock undertook to make a complete flexibility analysis 
of this piping configuration including all of the four secondary 
effects. The 1964 paper reports his results. It also reports that 
a computer program had been written to investigate this analysis 
and the resulting system of equations. The program had been written 
but was not developed sufficiently to thoroughly investigate the 
relationships involved. 

^ John E. Brock, "Some Secondary Effects In A Simple Piping Structure 
Under Heating," Journal of Applied Mechanics , Vol 31, March 1964, 
pp. 88-90. 
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1-2 Objectives of Work 

The object of this thesis was to write a computer program for 
a large digital computer that would be capable of performing the 
calculations required by Professor Brock's theoretical development 
and one that would allow application of the theory to any engineering 
problem that involves a piping structure of the given configuration. 
Further, it attempts to determine under what conditions these 
secondary effects may become significant in calculation. 

1-3 Assumptions and Limitations 

The investigations that are reported herein are based on an ell- 
shaped piping configuration with a square corner. No other shape is 
considered. Also, the configuration is assumed to have deformation 
only in the plane of the legs a and b. 

It is further assumed that deformations and displacements are 
small compared to the longitudinal dimensions of the pipe, except 
that a first approximation is included for the difference between 
arc length and chord length of the deformed pipes. 
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1-4 Notation 



a ,b 
A 




e 



E 



F 1' F 2 



g ij 



H 



I 

L 



M,N 



P 

P 

Q 

U 

V* 

vJ 



X 

X 

y 

Y 



TABLE I 
NOMENCLATURE 

length of straight pipe elements 
cross sectional area of pipe material 
cross sectional area of pipe contents 
unit thermal strain 
Young f s modulus of elasticity 
structural forces in ell-configuration 
coefficients defined in equations (16)- (19) 
shearing modulus of elasticity 

= pressure coefficient 
moment of inertia of pipe cross section 
length of pipe 
bending moments 
internal pressure 
lateral force 
axial compressive force 
Vj' /l 

(-ur)* 

ql 2 /ei 

axial coordinate 

horizontal deflection of joint 

lateral coordinate 

vertical deflection of joint 

coefficients (see Table II) 
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TABLE I (cont.) 



M* 

X 

8 

A 

7 

a 

A 

f 

f 

f 



coefficients defined in equations (15A) 



and (15B) 

CP/A6 = effective shearing strain 

tip axial extension 

tip lateral deflection 

conveniently small positive numbers 

shear distribution factor (see equation (9)) 

^ei/agl 2 



rotation of joint 

difference between arc and chord (first 
approximation) 

coefficients in equations (9), (10), and (11) 
Poisson's ratio 

(inside diameter of pipe) / (outside diameter 
of pipe) 
tip rotation 
xu/L 
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FORMULAS FOR (X COEFFICIENTS 
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SECTION II 



THEORETICAL ANALYSIS AND TECHNIQUE FOR INVESTIGATION 

2-1 Review of Flexibility Analysis for the Ell-Shaped Pipe 

Before presenting the investigation of this ell-shaped piping 
configuration in detail, it is felt that certain principles re- 
garding the analysis merit some review. 

If the tip loaded 

cantilever beam shown in M 

Figure 1 is analyzed by 
writing a moment equation 
about point A the following 
equation results: 

Figure 1 

< l > Ely" = M + Q(A-y) + PCl-x) 




Solving this differential equation and applying boundary 



conditions yields 

(2) ys |tApa(\-cos|)seca+ PlX ' ^ C i- Co$^)tanu] • WfJ' 

^ EX u? 



for boundary conditions 



a. y(0) = 0 

b. y(L) = A 

c. y ' (0) = Y 

d. y ' (L) = Cp 



See Table I for nomenclature. 



In addition 

o) a = CMd, + pl(c < 4 + ’)oIj)J(lVei') 
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(4) 



and 



cp * [Mc* 3 +PL(^,+ *)c< 4 )](l/ei) 



The alpha coefficients appearing in the above and subsequent 
equations are defined in Table II with alternate forms to facilitate 
computation. 

In terms of all four secondary effects, <5 is given by equation (5). 

(5) 5 = el - zjh >J= S % - - ?\ 

The first term on the right hand side of equation (5) represents 
thermal expansion, the second represents axial contraction resulting 
from hoop tension caused by internal pressure, the third represents 
axial contraction due to axial compression in the metal, and the 
last represents the effective axial contraction resulting from the 
difference between arc and chord. The pressure terms may be combined 
to yield the equation, 

(6) & = (e+f>M)L - Qh/AE - ?\ 



where ^ is given by ^ 

<» A = i j ((af-^)dx 

when higher order terms are neglected (see Appendix G) . If alpha 
coefficients are substituted, 



(8) A - (M l <x 8 + 2.0 MPLotj + P 2 eoQ( 
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Other pertinent relationships include the shear distribution 

^ 2 
factor 

( 9 ) £ = /*! + Alf'/Ct-hf'? 

where 

(10) /* i ' -A + 8if)/<b(/+lA) Z 

and 

(id /<* = (I6+Z0i>t 8 jJ z )/z (i t-Vf 

To apply the foregoing analysis to an ell-shaped pipe, the 
corresponding quantities are defined as they appear in that 
configuration. The basic structure is that shown in Figure 2. If 
the configuration is subjected to uniform heating, thermal expansion 
creates a situation resembling that of Figure 3. 





The resulting free body diagrams for the two legs a and b is depicted 
in Figure 4. 



John E. Brock, "Shear Distribution in Piping," Heating Piping and 
Air Conditioning . Vol 35, January 1963, pp. 141-143. 
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Referring to the horizontal leg as part 1 and to the vertical leg as 
part 2, Table III list the quantities corresponding to the general 
case for each part. 



TABLE III 



ITEM PART 1 PART 2 

Lab 



P 

Q 

M 




1 

-N 



A 

6 

f 



X 

Y 

9 



Geometrical continuity requires equating corresponding expressions 
for X,Y, and 9; i.e., 



6a. = A b 
6b = A<x 
9 *- cp b 
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(i2> ea- f ^ e - fi + M =[-^,4 l\b(^ + ^ 3 )]t>y EI 
d3) eb- ?i =[j^^i+F l a(o( 2L +») o( i)J a /£i 

(i4) A/aoC^- ^^(oC, +- 5^)' - Mb 5(^4 + 0^0 



A superior bar indicates that the quantity has been evaluated 
using the parameters of part 2 and no superior bar indicates evaluations 
with the parameters of part 1. 

Solving equation (14) for the moment, N, yields 
(15) 



Nc F ' a ' 1 '( < V9 0 W) (St, +>) 5^ 

ao^ + 



= F , + ^ a. 



where (15A) ^ + bo(j) 

and ( 15B) $t.s ■f bo^ 

Now define; 

(16) 3.1= -^,^^/eI 

07) 3**[b’^+H)^ 

(18) «j al r [ce(« t + *)oQ)+ hi - <*, ^ o>]/ EI 

(19) = 0<» aVEi 



and the system of equations become 



' 


3m 3 * 








' (e+p(») x - A”\ 


(20) 


J3». 3»\ 






mm 


_(e+*>H)b -7i_ 
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In order to solve for the inknowns F^ and F^ , an iterative 
procedure is strongly suggested by the nature of the equations. 

By assuming initial values of F^ and F^, which can be any 
value, one can calculate finite values for all the coefficients 
g _ and the quantities on the right of equation (20) . This equation 
is solved for new values of and F^ which in turn provide improved 
values for the next iteration. This procedure of iteration is 
continued until satisfactory convergence is reached. Once convergence 
is obtained such quantities as moment, 0, and other desired values 
are obtained using the converged values of F^ and F^* 

2-2 Dimensionless Parameters 

The large number of variables involved in this problem could 
produce millions of different problems. In order to pursue the 
investigation as thoroughly as possible and at the same time retain 
some compactness in the results, four dimensionless parameters were 
used. These are: 

(1) Leg Ratio a/b 

(2) (Radius of Gyration) / (Total Length) 

(3) Pressure /Modulus of Elasticity 

(4) Thermal Strain 

To cover reasonable ranges of all four parameters evidently 
requires solution of a large number of problems each of which is 
itself rather demanding. Accordingly, use of a large digital computer 
is definitely required. 

Each of the four secondary effects can either be included or 

Zj. 

excluded in the calculation. This produces 2 or 16 different possible 
combinations of secondary effects. 
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2-3 Requirements Imposed on Digital Computer Program 



A computer program for studying this problem has certain basic 
requirements . 

A. First of all, it must be capable of accepting certain input 
data, which includes a basic physical description of the 
problem as well as control variables that instruct the 
program in how to treat these data. 

B. From the input information it must recognize the parameter 
to be investigated and must load the various values of the 
parameter over the desired range. 

C. For each parameter value in the given range it must: 

1. Solve for the F's by the iteration procedure and then 
determine other desired values. 

2. Accomplish C.l. for all different combinations of the 
four secondary effects. 

D. It must increment the parameter value and repeat step C 
until results for all values in the range have been found. 

E. It must provide results in either tabular or graphical form. 

The following section describes the computer program in more 

detail and shows how the above listed demands are met. 
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SECTION III 



DESCRIPTION OF COMPUTER PROGRAM 



3-1 General Remarks 

As mentioned in Section II and suggested in reference 2, a 
computational procedure is evident. Assume initial values of 

and F^ (which may be taken to be zero) , evaluate the coefficients 
g_ and the terms and ^ and upon salving the system of equations (20), 
obtain improved values of F^ and F^, continuing until satisfactory 
convergence is reached. This idea along with the other basic require- 
ments (A-E of Section II) composes the foundation upon which the 
computer program is constructed. 

The program was first written in Fortran 60 computer language 
for use on the Control Data Corporation Model 1604 digital computer 
located at the Naval Postgraduate School, Monterey, California. In 
the midst of the investigation reported in this paper, the CDC 
computer was replaced at the Postgraduate School by an International 
Business Machine Model 360/67. The present version of the program 
is written in Fortran IV language for use on the latter machine. 

All computations in the program are executed using the Double 
Precision Capability of the system. This feature improves the 
accuracy of the results but unfortunately requires more time than 
the single precision calculation. 

An effort was made to keep the program as general and as flexible 
as possible. All of the information needed to specify a given problem 
is read directly from data cards. The information on the data cards 
defines the physical condition to be investigated and gives the program 
additional control information concerning the parameter to be used 
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in evaluation, range of the parameter, combinations of the secondary 
effects to be investigated, convergence criteria, etc. 

It allows the user four different parameters with which to 
evaluate the problem. It allows for all sixteen possible combinations 
of the four secondary effects and computes all sixteen of these for 
each value of the parameter. 

It has optional output. In the data cards the user specifies 
the desired type of output: tabular, graphic, or both. 

To describe the computer program more completely requires 
explaining in more detail the functions of its seven major parts. 

A. Input 

B. Loading of Parameter Values 

C. Consideration of All 16 Different Combinations 

D. Computation by Iteration 

E. Incrementing Parameter to Next Value in Range 

F. Output 

G. Subroutines 
3-2 Input 

The term input is used here in referring to the information that 
is supplied to the program via the data cards. In all, there are 
fifty-five different items read in as data. These items completely 
define the physical situation under consideration. Further, it informs 
the program of the dimensionless parameter to use in the investigation, 
the range of the parameter, conditions to investigate, and so on. 

Each of the variables read in as data is listed below along with a 
description of its function in the program. The order of this listing 
is the order in which they appear on the data cards. The names of the 
variables are the names which are used as input. The same physical 
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quantity may be called by other names later in the program and it 
should be remembered that the variable name mentioned below is 



applicable only in the input section of the program. 



R 


length of the shorter leg, a, (inches) 


S 


length of the longer leg, b, (inches) 


ZIN 


pipe inside diameter (inches) 


ZOUT 


pipe outside diameter (inches) 


DTEMP 


temperature difference (degrees fahrenheit) = 
final initial 7 


Q 


internal pressure (psi) 


E 


modulus of elasticity of pipe material (psi) 


COE FI 


coefficient of linear expansion of pipe material 
(in/in - °F) 


KPARl 


Integer variable which can assume the integer values 
1 to 4 inclusive, that designates the parameter 
to be used in the investigation as follows: 

KPARl = 1 = Leg Ratio 

KPARl = 2 = (Radius of Gyration) /(Total Length) 

KPARl = 3 = Pressure/Modulus of Elasticity 
KPARl = 4 = Unit Thermal Strain 


L 


Integer variable that specifies the size of the major 
parameter increment (see Section 3-3). This can be 
any value but for graphic output should not require 
more than thirty stored values of the parameter. 
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Integer variable that designates the initial value 
of N, itself an integer variable which is used 
as a subscript for the parameter values. This 
(the value of NN) is useful if the program ever 
exceeds the time allotted for one run and it is 
terminated by the operator, or if the given size 
of increments do not yield convergence over the 
desired range. By giving NN the value of the last 
value of N on the terminated run, the user is able 
to modify certain other control variables such as 
L, LL1, FTE's, etc., for the next computer run and 
begin computation at the point where the last run 
ended, regardless of the reason for termination. 

For the first run of a given parameter or problem 
the value of NN is always 1. 

Integer variable that designates the final value of the 
subscript N in the parameter range. 

Integer variable which designates the initial value of 
MT. The significance of MT is explained in more 
detail in the section defining loading. MB can take 
on any value from 1 to 16. 

Integer variable that indicates the final value of MT. 

MC 2: MB always and 1^MC^16. Both of the above 
conditions must be satisfied. The program is designed 
in such a way that if MT does not equal its final 
value, MC, then MT will be incremented by 1 , or the 
next combination of secondary effects will be 
calculated using the same parameter value. What 



this in effect does is have the computer make all the 
calculations for all the combinations between and 
including MT = MB and MT = MC. For example, if 
MB = 3 and M.C = 7 for a given set of input values, 

MT would take on values 3, 4, 5, 6, and 7 for each 
value of the parameter. Consequently, the output 
would reveal results for the combination correspond- 
ing to the 5 different values of MT. 

MD Integer variable that designates the type of output 

desired . 



1 - 


MD 3 




MD 


= 1 


Tabular output only 


MD 


= 2 


Graphic output only 


MD 


= 3 


Both tabular and graphic output 



MIT Integer variable that indicates the limit to the 

number of iterations the program will perform. If 
this number is reached before the convergence 
criterion is met, the program will terminate 
operation on its current set of data, give the 
necessary output, then proceed to the next set 
of data. 

KKK Integer variable that assigns a value to the constant 

IA. IA is used in a relationship that determines 
the improved values of F's to be used in the next 
iteration on the basis of the last assumed values 
of the F's. The relationship involved is: 

Fl-= 4 » Fl(Hew) IA* PI(qid) 

4 + IA 
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and a similar relationship exists for the improved 
F2. It was found that generally as the value of N 



LL1 

EP2 



EXSCAL 

Y SCALE 
IXUP 

IYRIGH 

FTE11-FTE 162 



increases IA should increase also, thereby making 
the improved value closer to the old value than 
the new value. 

Integer variable that specifies the size of the sub- 
increments of the parameter. (See section concerning 
loading) . 

This value sets the convergence criteria. After each 
iteration the new values of the F's are compared to 



- Fl(oLO) 


L±! 




2(010)1 


\) (Fi(WEW))* +- ( 


F2/N£vV)P 








If this is satisfied the values are considered converged. 

This quantity is the argument used in the DRAW subroutine 
(see Appendix A-4) that assigns the scale to the X-axis 
or abscissa axis when graphic output is desired. 

Same a EXSCAL above except it governs scale on Y-axis. 

Another argument for the DRAW subroutine that specifies 
how many inches up the graph one wants the X-axis 
located . 

Specifies how many inches from the left side of the graph 
one wants the Y-axis. 

The 32 different values are the last converged values 
of FI and F2 for a particular combination of secondary 
effects. The final digit distinguishes between FI 
and F2. The other numerical digits run from 1 through 
16 and correspond to values of MT. These FTE values 
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are retained by the computer and used as the first 
approximation for the F's in the corresponding 
combination when the parameter has been incremented. 

In most investigations the values assigned by the 
data card are all zero. These values, however, are 
just the initial values assigned to the FTE's. Once 
computation begins the computer takes over the task 
of assigning these values. Every time convergence is 
reached the values of two FTE's change and as 
convergence is reached for all combinations and one 
parameter value, all 32 change their values. Use of the 
variables aids the computation process in two ways. 

The nature of the governing equations causes difficulty 
with convergence under some conditions. Therefore, 
if the first assumed values of the F's are near the 
correct values a better chance for convergence exists. 
If, on the other hand, the guess is arbitrary, or not 
close to the converged value, the system of equations 
may not converge. If the last converged values are 
assumed as the first approximations with the incremented 
parameter, they should be good estimates depending on 
the size of the parameter increment. The smaller the 
parameter increment the better the estimate should be. 

Another utilization of the 32 variables allows 
the user the capability of picking up computation at 
any point in a parameter range. Suppose, due to excess 
time the run was terminated before the entire range 
was completed. Then, on the subsequent run, the user 
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would be able to start at the same place in the parameter 
range that the previous run ended. This provides 
continuity without duplication and the ability to read 
reasonable first approximations to the F's in as data 
regardless of where in the range one desires to continue. 

3-3 Loading of Parameter Values 

The integer variable KPAR1, read from the data card, instructs 
the program which dimensionless parameter will be used to investigate 
the structure. For a given set of data (one set of six data cards) 
only one parameter can be investigated. Each of the four values that 
KPAR1 may take on indicates a different loading scheme for the six 
subscripted variables. Once the program knows the value of KPAR1 it 
proceeds to the proper loading scheme for that parameter. 

The six subscripted variables involved are: 



AA(N) 


leg a 


BB(N) 


leg b 


XDIN(N) 


pipe inside diameter 


XDOUT(N) 


pipe outside diameter 


XP (N) 


i nternal pressure 


ALIT(N) 


unit thermal strain 



Here again the variable names associated with the physical 
quantities apply only in the loading scheme. The same physical quantity 
may have another variable name at another place in the program. See 
Appendix A for a complete listing of computer program nomenclature. 

The subscript N takes on integer values from 1 to MA. For each 
value of N, the six subscripted variables have values stored in the 
computer memory. In each of the different loading schemes only one 
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of the six subscripted variables takes on different values for each 
different value of N. The other five have fixed values. 

Associated with N in the loading procedure are two other integer 
variables, MA and L. The values of these two quantities determine 
the range of the parameter and the size of the main parameter increment. 

An important feature of the loading process is that the parameter 
values are loaded in such a way as to produce the least critical loading 
conditions on the structure for the parameter value associated with N = 1 to 
the most critical loading condition when N = MA. This essentially, 
means that at the least critical end of the range the initial assumption 
of the arbitrary values of the F's will probably yield convergence. How- 
ever, as the parameter proceeds to the more critical end of the range, 
the facility utilizing the FTE's becomes more and more necessary in 
achieving convergence . 

3-4 Combination of Secondary Effects 

The capability that allows the computer program to solve the system 
of equations for all sixteen different combinations of secondary effects 
for each parameter value is afforded through the subroutine FACTOR. 

To supply the subroutine with information it needs, the integer variable 
MT is used. MT takes on values 1-16 for each parameter value in the 
range. Each of the 16 values represents a different combination of 
secondary effects. When the subroutine is called with MT specified, 
it responds by supplying the main program with four factors and a 
particular designator. These four factors are multiplied by the terms 
in the system of equations that represent the four different secondary 
effects. Each of the factors can have only two different values: 1.0 
and 0.0. It is therefore evident that if the combination requires the 
suppression of an effect in the calculation, the governing factor is 
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zero; on the other hand if the effect is to be included the 
governing factor is unity. 

The main program is so designed that all of the sixteen different 
combinations do not have to be calculated each time, say in a situation 
where the user is only interested in one or two combinations. If such 
is the case the integer variables MB and MC are manipulated to provide 
the combinations desired. MB sets the initial value of MT. After 
each combination is completely solved for all desired quantities the 
main program tests the value of MT against that of MC . MC represents 
the final value of MT to be considered for each parameter value and 
in the test if MT is less than MC its value is incremented by 1 and 
that combination is investigated. This process continues until 
MT = MC and after the MT = MC combination is solved the program increments 
the parameter. The limitation on this capability lies in the fact that 
combinations represented by values of MT between MT = MB and MT = MC 
will be calculated even if some of these intermediate cases are not 
specifically wanted. 

The designator referred to earlier is realized as a label in the 

printout of the results. It is printed in a certain column of the 

results and is a four character identifier that tells the value of 

each factor in that particular combination (see the explanation of this 

subroutine for specific meanings for these designators in Section 3.8). 

3-5 Computation by Iteration 

The iterative procedure, as mentioned earlier, is that of assuming 
initial values of FI and F2 (which may be taken to be zero), evaluating 
the coefficients g_ and the terms ^ and and upon solving the system 
of equations (20), obtaining improved values for the F's. This is 
continued until satisfactory convergence is reached. Through the use 
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of the values of the FTE's an attempt is made to minimize the number 
of iterations required. This saves computer time as well as increases 
the probability of convergence. 

After each iteration the newly calculated values of the F's are 
tested against the old ones using the following convergence criterion: 



|fi(hew) - FI (olp)| -(• [ (NEW) -tZ (0<-O)| 

V(RU*'ew)) i + (f 2 (ne w)f 



£ EP Z 



If this criterion is not satisfied, the program selects improved 
values for the next iteration and sends these newly selected assumptions 
of FI and F2 through the system of equations again. The process for 
the selection of the new F's is defined by the following two relation- 
ships : 

FlNPXT g 4* FlN + JA* FI . F2.NEKT s 4-xFZN +IA* 

4 4-IA > 4 + XA 



where F1NEXT and F2NEXT are the improved values of FI and F2 respectively 
to be used in the next iteration. FIN and F2N are the new values of 
F^ and F^ just computed and FI and F2 in the equations are the values 
of F^ and F^ assumed at the beginning of the iteration. The factor 
IA is an integer whose value is set by the integer variable, KKK, 
read in from the data card. It is obvious that the improved value can 
be closer to the older value or the newer value at the discretion of 
the user. The nature of the equations is such that when the parameter 
value produces less critical loading conditions, convergence is achieved 
more rapidly if the improved value is taken to be closer the new values. 
Conversely, when the parameter is in a more critical region, a selec- 
tion of an improved value closer to the old value produces convergence 
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more rapidly. The term ’’critical region” refers to that range of the 
parameter where the system of equations behave in a highly nonlinear 
fashion. 

There is no set rule for this selection of improved values and 
the best instruction is gained through experience with the program. 

Using values of IA between 0 and 10 is normally adequate for most 
situations. 

If the convergence criterion is satisfied, the latest calculated 
values of the F’s are accepted as the correct ones for the given value 
of the parameter and the given combination of secondary effects. 

Using these values of the F's other necessary quantities are calculated. 
The program then goes to the next combination of secondary effects 
and begins iteration for the solution of the F’s for that situation. 

The first approximations for the F’s are the converged values found 
for that combination with the previous parameter value. When the 
program has solved the equations for all combinations desired, it 
increments the parameter and the process starts all over again. 

3-6 Incrementing the Parameter 

Two different parameter increments are utilized in the program. 

They are referred to hereafter as &PAR and <§par. The two can be 
described as follows: 

If PAR(N) is the parameter for a particular value of the subscript, 
then PAR(N + 1) = PAR(N) + &PAR. Therefore, Z^PAR is the increment 
whose size is the difference between two consecutive stored parameter 
values. The size of &PAR is determined by appropriate input data 
and by the control variable L (see description of L in Section 3-2). 
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The sub- increment , £par is used to give the program the capability 
of making the parameter as small as deemed necessary in order to achieve 
convergence. The basic theme of the investigation of any of the four 
dimensionless parameters is to proceed from that end of the range 
that provides the least critical load toward the end which produces 
the most critical condition in the range. Each time the parameter 
is incremented, the last converged values of the F's are used as 
the first approximation. By making the increments smaller, the first 
approximation is closer than if the parameter had a larger increment. 
Therefore, it is conceivable that by making these increments small 
enough, convergence can always be realized. It is for this reason 
that the program uses increments of two different sizes. 

For the graphic output capability, it is necessary to store many 
values of the variables that are used in the graphs. The storage 
capacity of the computer thus establishes a lower limit for the larger 
increment. In order to further decrease the size of the increment the 
smaller increment, <$par, is used. Except when they coincide with 
APAR values, parameter values resulting from this smaller increment 
and the corresponding calculated values are not stored, but used only 
to obtain new FTE's for the subsequent parameter value. Another way 
of expressing the process is to say that after all necessary calculations 
are made for one parameter, the parameter value is incremented by £par. 
Using the last converged values of the F's as first approximations, 
converged values for the new parameter value are sought. Thus the 
incrementation continues over the parameter range. If the first 
parameter value considered was a stored value, then after LL1 sub- 
increments the next stored parameter value will be reached. LL1 is 
the integer variable read from the data card. 
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Figure 5 aids in understanding the use of the two increments. 



^(P&r&me+er) 




J I ! ! l_ 

PARAMETER. 



Figure 5 



From Figure 5 the following relationships are apparent: 

LL1 x & par = APAR PAR(N + 1) = PAR(N) + APAR 

Also, for the parameter values corresponding to the solid vertical 
lines, necessary calculated values are stored for graphic or tabular 
output. For parameter values corresponding to the broken vertical 
lines only converged values of the F f s are calculated for use as 
the first approximation when the parameter is incremented. 

To summarize, each time the parameter is incremented it is 
increased by the amount £par. Every LL1 times, all calculated values 
necessary for output are stored. The size of APAR is somewhat 
restricted by the storage capacity of the computer but the size of 
(5 par has no restrictions. It must be remembered, however, that by 
decreasing the size of £par, the computing time for a given range of 
the parameter is increased accordingly. Therefore, in the interest 
of time on the computer, the value of £par should be as large as possible, 
commensurate with providing the necessary convergence. 

3-7 Output 

The user has some choice concerning the type of output. The 
results can be represented in tabular form, graphical form, or both. 
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To specify the type of output desired for a given set of data^ assign 

the proper value to the input integer variable MD. MD can take 

on values 1 to 3 with the following significance attached to each: 

MD = 1 Tabular Output only 

MD = 2 Graphic Output only 

MD = 3 Both Graphic and Tabular Output. 

If MD is set equal to 1, ten columns of calculated values will 
be printed out for all of the desired combinations of secondary 
effects. The values shown are: the value of the subscript (N) , 

the value of the integer MT, the combination designator, moment, 
moment ratio, the converged values of F^ and , the value of the 
parameter, and the number of iterations required for convergence 
in each case. In addition at the top of each page of output are 
listed the values of the raw data exactly as it is read in from the 
data card. These values only serve as a convenience in showing what 
was read Into the program and what the individual values were prior 
to the parameter loading. Other values that are printed include the 
exact moment, the parameter that was investigated, the elapsed time 
in seconds since the last stored parameter value with its associated 
values were calculated, and the critical loads for each parameter 
except the leg ratio. Actually the output sheet is very straightforward 
if one follows the column labels. 

Some of the above mentioned terms deserve a little more explanation. 

Combination Designator: The four character designator supplied 

to the main program by the subroutine Factor. The four characters are 
either zero or one and they represent the value of a particular factor 
for the given combination of secondary effects. For instance, if the 
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designator is 1001, the factor that is multiplied by the term represent- 
ing the secondary effect, arc and chord difference, is 1. The factor 
that is multiplied by the secondary effect term for shear is 0. The 
factor multiplied by the term representing axial deformation is 0 and 
the factor multiplying the beam column term is 1. The order of the 
factors always appear in the same order. Arc and chord, shear, axial 
deformation, beam column action appear in that order from left to right. 

Exact Moment: When the term exact is used referring to a calculated 

quantity in the system of equations it means that the calculation was 
made with all four of the secondary effects allowed and none of the 
effects suppressed; for such a calculation, the combination designator 
is 1111. 

Moment Ratio: Ratio of moment calculated with given combination 

designator to exact moment. 

If MD is set equal to 2, only graphic output results. For graphic 
output certain limitations exist: 

a. The ordinate is always the moment ratio. 

b. The abscissa is always the parameter under investigation. 

c. MA or the number of stored values for each subscripted 
variable should not exceed 30. 

One graph is obtained for each set of data investigated. On each 
graph six sets of values are plotted. These six sets correspond to the 
first six values of MT. See Subroutine Factor in Section 3.8 for the 
combinations of secondary effects that correspond to the six values. 

The combinations included in the six are with all effects allowed (exact 
solution), all effects suppressed, four combination where only 1 of 
the four effects is allowed and the other three are suppressed. There- 
fore, if the user specifies graphic output only, the input integer 
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variables MB and MC should equal to 1 and 6 respectively. Calculations 
for all other combinations would be redundant since there is no 
tabular output and the first six combinations are all that will be 
represented on the graph. If on the other hand, tabular output is 
selected or the user desires both types of output on the same run 
then the above restriction is not applicable. 

Before one employs the graph plot form of output it is necessary 
to have some knowledge of the values that are to be plotted. The 
data card contains four of the arguments to the DRAW subroutine. 

For a meaningful graphical output a reasonable scale has to be 
selected and other parameters have to be chosen. The recommendation 
therefore is to use the tabular output first and from it plan the 
graphical output so as to describe the results most usefully. In 
addition the user should be familiar with the information relating to 
the DRAW subroutine contained in the next section and in the subroutine 
system library write-up reproduced in Appendix A-4. 

3-8 Subroutines 

There are three subroutines used in the program. The first two 
subroutines, designated ALPHA and FACTOR were programmed by the author 
and the third, designated DRAW, is from the subroutine library at 
the computer facility of the Naval Postgraduate School, Monterey, 
California. The complete library write-up forDRAW and its use appears 
in Appendix A-4. 

ALPHA: This subroutine is used by the main program to compute 
the values of the alpha coefficients as they are shown in Table II, 
page 13 . It is called twice during each iteration and supplies the 
program with values of the corresponding to the information 

that is given to it, namely, , >0 , and EP1 or Tjl 1,r\ , and EP1 , 
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depending on whether the main program wants the values of 

or 0( S . 
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FACTOR: This subroutine is used by the main program to set 

the combination of secondary effects to be considered. As previously 
mentioned, it works along with the integer variable MT which takes 
on values 1 through 16. Each of the values designates a certain 
combination of the secondary effects. Upon receiving this number 
designation, the subroutine supplies the program with four factors 
each of which has the value of either 1.0 or 0.0 depending on the 
combination. These factors are multiplied by the quantities that 
introduce the secondary effects into the system of equations. If 
the factor is 1.0 that particular secondary effect is included in 
the calculations; if it is 0.0 the secondary effect is suppressed in 
the calculation. Tables IV and V on the following two pages give more 
information concerning these combinations and their designations 
that is helpful in understanding the particular system for varying 
the combinations. 

DRAW: This subroutine is the one through which the program 

affords the user the graphic output option. The complete write-up 
of this subroutine just as it appears in the Postgraduate School’s 
subroutine library is reproduced in Appendix A-4. It is felt that 

no further explanation is necessary except for the four arguments 

0 

that are read in as data in the main program. Referring to the 
write-up in Appendix A-4, the arguments h. through k. are supplied 
to the calling statements by the last four values on the second 
data card. This allows user control of the scales desired and the 
axis location. 



38 



TABLE IV 



FACTOR NAME 


SECONDARY EFFECT THAT IT CONTROLS 


FAC1 


Difference in Arc and Chord 


FAC2 


Shear 


FAC 3 


Axial Deformation 


FAC4 


Beam Column Action 
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DESIG- STATUS OF EACH SECONDARY EFFECT VALUE OF EACH FACTOR 

NATOR 

A-C SH AX-DEF. B-C FAC1 FAC 2 FAC 3 FAC 4 
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A means effect is allowed 
B means effect is suppressed 
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SECTION IV 



RESULTS OF INVESTIGATION 



4-1 Errors in Paper No. 63--APMW-21 

Prior to developing the computer program to utilize the relation- 
ships of Section II, certain review and checking was undertaken by 
the author. This review was undertaken principally to gain a working 
familiarity with the equations and nomenclature. During the course 
of this review and later in developing the program a few errors 
were detected in the original paper. 

Most of these mistakes were algebraic in nature. One was in 
the area of the critical buckling loads for beam column considerations. 
The mistakes are: 

1. Equation (8), the factor "2.0" in the second term of the 
numerator did not appear in the same equation in Brock’s paper. 



There were three mistakes in the formulas for the o( coefficients 
of Table II. 



5. The last incorrect item found was in the statement relating 
to the critical buckling loads. The paper indicates that 
defines the critical buckling load for the ell-shaped pipe structure. 
This is true only for the equal leg case but is not true in general. 

All of the results reported in this paper reflect the corrections 
as stated above. 



2. & 3. In the formula for Ot ^ , the 0^5 in the first two terms 
of the series should be In the reference paper they are O^S . 



4. The formula for , the last term should 
opposed to " in the paper. This was apparently a typographical 

error. 




41 



4-2 Numerical Calculations 



With the use of the digital computer program described in 
Section III, calculations were made for several sizes of pipe 
and for various values of temperature increase, as well as for 
different values of the other inputs. This was done not only to 
test the integrity of the computer program, but also to obtain 
sufficient data to permit a search for trends and useful modes 
of presentation. It turns out that the leg ratio, a/b, is the 
most significant parameter, as might have been expected, and 
presentations are confined to displaying the moment ratio as a 
function of leg ratio for various combinations of secondary 
effects all curves shown being for the following specific case. 



Length of leg a 
Length of leg b 
pipe inside diameter 
pipe outside diameter 
Temperature difference 

pipe internal pressure 

pipe modulus of Elasticity 

pipe coefficient of linear 
expansion 



varied 

360 inches 

10.020 inches 

10.750 inches 

530 degrees fahrenheit 
(70-600 degrees F) 

400 psi 

27.4 x 10^ psi 

7.23 x 10 6 in/in deg F 



4-3 Final Results and Conclusions 

The results are shown in the form of eight curve sheets, 

Figures 6 through 13, which should be considered in pairs, the first 
of each pair showing the effect of suppressing (i.e, failing to 
consider) only a single secondary effect and the second showing the 
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effect of suppressing the other secondary effects and considering 
only this one. For comparison purposes, each curve sheet also 
shows the effect of neglecting all four secondary effects. 

These curves are based only upon what has been called the 
ll moment ratio 11 , i.e., the ratio of the moment at the intersection 
of the two legs as calculated by an inferior procedure to that 
calculated by the best available procedure, namely that of including 
all four secondary effects. It is clear that substantially the same 
sort of result would have been obtained by considering a similar 
moment ratio for either of the two built-in ends or a similar ratio 
of reaction force. 

These curves may be regarded as speaking for themselves. The 
most significant conclusion that may be made is that neglecting 
secondary effects, indeed can have a profound influence on the 
accuracy of piping flexibility calculations, leading to gross errors 
in the evaluation of stresses in cases of leg ratios which, though 
extreme, are likely to be encountered in piping practice. 

Since the secondary effects are quite subtle to evaluate, this 
conclusion could not have been forthcoming from a less sophisticated 
and demanding approach. It is believed to be an entirely new con^ 
elusion and it is certainly a significant one. 
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APPENDIX A 



COMPUTER PROGRAM NOMENCLATURE 



A-l MAIN PROGRAM 

NOTE: Where possible the variable names of the computer program 

are identified by the notation of Tables I and II. Also, in some 
cases the same physical quantity may be represented in the program 
by more than one variable name. In such cases the extreme right 
hand column will indicate the general area of the program where the 
physical quantity takes on the indicated variable name. 



VARIABLE NAME 


TABLE I & II NOTATION 


COMMENT 


A 


a 


Computation Section 


AA(N) 


a 


Loading Section 


A1-A10 


Of,— Of,o 




AB1-AB10 


Of, — Of/o 




ALIT(N) 


e 


Loading Section 


ALPHA 


Subroutine ALPHA 




AREA 


A 




B 


b 


Computet ion Section 


BB (N) 


b 


Loading Section 


B1 


fetpH)a.- A 




B2 






BET1 


L®- 1 fa , + l ) o< +)T/ (<*•&} + 




BET2 


C b 1 (of , + 9 » • + b a, 1 ) 


COEFI 


Material coefficient of thermal expansion 


DELT1-DELT4 


Size of 6par for each parameter 1-4 respec 


DET 


3«i * 9 m. “ x ^ 


III 


DIN 


pipe inside diameter 


Computation Section 


DOUT 


pipe outside diameter 


Computation Section 
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COMPUTER PROGRAM NOMENCLATURE (cont.) 



VARIABLE NAME 


TABLE I & II NOTATION COMMENT 


DTEMP 


Difference of Temperature 


E 


E 


ELIT 


e Computation Section 


EP1 ,EP2 

ETA 

ETAB 


e 'A 

0 

9 


EXMOM(N) 


Exact moment for a particular parameter variable, i.e. 



with all secondary effects allowed. 



EXSCAL 


Applicable to graphic output: X scale in units per inch 


FAC1-FAC4 


Four factors used for either suppressing or allowing 



their corresponding secondary effects. 



FACTOR 


Subroutine FACTOR 


FI ,F2 


Assumed values of the forces F^jF^ at the beginning 



of an iteration. 



F 1CRIT , F2CRIT 


Euler Buckling load for F^ and F^. 


F1N,F2N 


Newly calculated values of F^ and F^ from a particular 



iteration. 



FTE11-FTE162 


See section on input Data 


G 


G Computation Section 


8 ij 


coefficients defined in equations (16)- (19) 


H 


H 


I 


Counter use to count the number of iterations required 



for convergence. 



IA 


Factor to aid in determining improved values of F's 


LABEL 


Argument for Subroutine DRAW 


ITITLE 


Argument for Subroutine DRAW 
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COMPUTER PROGRAM NOMENCLATURE (cont.) 



VARIABLE 

IXUP 

IYRIGH 

J 

KHAT 

KKK 

KPAR1 

L 

LABEL 

LL1 

MA 

MB 

MC 

MD 

MIT 

MT 

N 

NN 

P 

PAR2 



NAME TABLE I & II NOTATION COMMENT 

Argument for Subroutine DRAW 
Argument for Subroutine DRAW 

Counter used to count the number of sub increments between 
the larger main increment. 

Secondary Effect combination designator 

Input variable that assigns value to IA 

Input variable that designates parameter to be investigated 

Input variable that aids in determining ^PAR. 

Argument for Subroutine DRAW 

Input variable that sets size of sub increments , £par. 

Input variable that designates final value of subscript, N. 

Initial value of MT 

Final value of MT 

Input variable that designates type of output 

Input variable that specifies the maximum number of iterations 
that will be performed before terminating operation on a 
particular set of data. 

Integer variable that specifies the combination of secondary 
effects to be considered. 

Subscript, used to identify calculated values of desired 
quantities with the corresponding parameter value. 

Initial value of the subscript N 

p Computation Section 

A designation for the second parameter, (Radius of 
Gyration) / (Total Length), in the computational part 
of the program. 
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COMPUTER PROGRAM NOMENCLATURE (cont.) 



VARIABLE 


NAME TABLE I & II NOTATION 


COMMENT 


PIE 


IT 




PIESQ 


TT l 




POIS 


y 


Computation Section 


Q 


P 


Input Data Section 


R 


a 


Input Data Section 


RAGY 


Radius of Gyration 




RHO 


f 




s 


b 


Input Data Section 


SFAC 


? 




w 


u 




WBAR 


u 




XDIN(N) 


pipe inside diameter 


Loading Section 


XDOUT(N) 


pipe outside diameter 


Loading Section 


XLA 


a 




XLAB 


A 




XMU1 


H> 




XMU2 


Hi 




XN 


Moment in Computation 


Section 


XP (N) 


P 


Loading Section 


Y SCALE 


Applicable to graphic output: 


Y scale in units per inch. 


XI 


I 




ZIN 


pipe inside diameter 


Input Data Section 


ZLEGR(N) 


First parameter, a/b 




ZLEN 


Total Length, a + b 




ZMOM 


moment 


c 
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COMPUTER PROGRAM NOMENCLATURE (cont.) 



VARIABLE NAME 


TABLE I & II NOTATION COMMENT 


ZOOM 


Initial thermal strain 


ZOUT 


pipe outside diameter 


ZPAR2 (N) 


Second parameter, (Radius of Gyration) /(Total Length) 



ZPAR3(N) Third parameter, Pressure/Modulus of Elasticity 
ZRAT1(N)-ZRAT16(N) Ratios of the Moments of the sixteen combinations 

divided by the exact moment for that combination. 
ZRATIO Same as the ZRAT()*S except this value is not stored. 

A- 2 SUBROUTINE ALPHA 



VARIABLE NAME 


TABLE I & II NOTATION COMMENT 


C 


9 or 9 


EP1 




T1-T10 

TX 

Y 

Z 


Q(/ — O (, 0 <*, — O { )0 

Iwl or lufl 
(-'Ujf* or (•'&)*' 
(lit)* or (Hs)^ 



A- 3 SUBROUTINE FACTOR 



VARIABLE NAME 


TABLE I & II NOTATION COMMENT 


JT 


Integer variable, MT 



LLAT Secondary effect combination designator 

R1-R4 Four factors used for either suppressing or allowing 

their corresponding secondary effects. 
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A-4 SUBROUTINE DRAW-DE SCRIPT ION FROM SUBROUTINE LIBRARY 



A. IDENTIFICATION: 



TITLE: 



General Graph Output Subroutine 



SUBROUTINE NAME: DRAW (J7-NP S- DRAW for CDC 1604 revised for IBM 360) 



B. PURPOSE: 

This subroutine generates data for plotting graphs on a CalComp 
plotter. Provision is made for curve drawing and point plotting, 
automatic scaling, graph titling and axis annotating. An attempt 
was made to provide a considerable amount of flexibility, at the 
expense, necessarily, of a relatively large number of arguments 
and a rather high memory requirement. 

C. USAGE: 

1. Definitions: 

In what follows the word "graph" will be taken to mean one 
piece or frame of graph paper on which there may be plotted 
one or more curves and/or sets of points. A "curve" will 
mean a continuous line generated by the sequence of straight 
lines joining successive points of the set defining the 
curve. A "point plot" will describe the representation of 
a succession of points by means of symbols (such as a cross) 
on the graph. The points are not connected in a point plot. 

2. Calling Arguments: 

All necessary information is transferred to DRAW through the 
calling arguments. The call statement is: CALL DRAW (NUMPTS, 



DATE: 



July, 1967 



PROGRAMMER: 



Revised by Patricia Johnson 



X , Y ,M0DCUR , I TYPE , LABEL , ITITLE , EXSCAL ,Y SCALE , IXUP , IYRIGH , 



MODXAX ,M0DYAX , IWIDE , IHIGH , IGRID ,LAST) 
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It is important to realize that one and only one curve or 
set of points is plotted each time DRAW is called. However, 
it is possible to call DRAW repeatedly if several curves 
and/or sets of points are wanted on one graph. 

The calling arguments are as follows: 



a . NUMPTS : 


The number of points defining a curve (2 £ NUMPTS £ 900) , 
of the number of points to be point plotted (2 - NUMPTSSL30) 


b. X: 


The array of X-ordinates, X, must be dimensioned at 
least equal to NUMPTS in the calling program. 


c. Y: 


The array of Y-ordinates, Y, must be dimensioned at 
least equal to NUMPTS in the calling program. 


d. MODCUR: 


Controls the number of curves, and/or sets of points 
on one graph: 


= 0 


This is the only curve, or set of points, to be 
plotted on this graph. 


= 1 


This is the first of two or more curves, and/or sets 
of points, to be plotted on this graph. 


= 2 


This is an intermediate curve, or set of points. 


= 3 


This is the last curve, or set of points, for this graph. 


e . ITYPE : 


Controls the type of plot (i.e., curve or point plot): 


= 0 


This set of points is to be represented by a curve. 


= 1 


These points are to be plotted with a cross (x) . 


= 2 


These points are to be plotted with a plus (+) . 


= 3 


These points are to be plotted with a square (□) . 


= 4 


These points are to be plotted with a diamond (O). 


= 5 


These points are to be plotted with a triangle (A). 
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f. LABEL: This is a Hollerith curve or point identifier. 

If a curve is being drawn, LABEL must have 4 char- 
acters (including any blanks), and these will be 
reproduced beside the end of the curve. This 
argument can be set in the calling program by a 
declaration statement, e.g., 





REAL IABEL/4H1234/ 




or 


REAL LABEL/ '1234 '/ 




or 


REAL LABEL/4H / 


(The latter must be 
used when no label is 
wanted. ) 



If a set of points is being plotted, LABEL is an 8- 
character identifier. The first 4 characters will be 
reproduced beside the first point, and the last four 
characters will be reproduced alongside the last 
point. This argument can be set by a declaration 
statement, e.g., 

REAL* 8 LABEL/ 8HFRSTLA ST/ 
or REAL LABEL ( 2 ) / 4HFRST , 4HLA ST / 

or REAL* 8 LABEL /8H / (The latter 

must be used 
when no label 
is wanted) . 

The above arguments, a. through f. (and q.), have 
meaning every time DRAW is called. On the other 
hand, the remaining arguments, g. through p. , have no 
meaning except when MODCUR = 0 or 1. 



60 



g. ITITLE : 


An array of 96 EBCDIC characters, the first 48 of 
which will form the first title line, and the last 
48 the second. The array must be dimensioned in 
the calling program, must contain the user’s iob 
identification, and must have unwanted characters 
set to blank. This may be done with a declaration 
statement such as: 

REAL* 8 ITITLE (12) / f 96EBCDIC character 1 / 


h. EXSCAL: 


X-scale in units per inch. EXSCAL will always be 
rounded off to one figure significance. If EXSCAL = 0, 
the X-scale will be computed by DRAW. This is called 
auto-scale . 


i. YSCALE : 


Y-scale in units per inch, with properties corresponding 
to those of EXSCAL. 


j . IXUP : 


Distance, in inches, of the X-axis from the bottom of 
the graph (O^IXUP ^IHIGH) . This argument will be 
ignored unless MODXAX = 2, see below. 


k. IYRIGH : 


Distance, in inches, of the Y-axis from the left of 
the graph (0 ^ IYRIGH £IW IDE) . This will be ignored 
unless MODYAX = 2, see below. 


1. MODXAX: 


Determines the mode of the X-axis location: 


= 0 


The X-axis will be located automatically by DRAW, 
with the origin of Y on the graph. 


= 1 


The X-axis will be automatically located by DRAW, 
with the origin of Y moved (in one’s imagination) an 
integer number of inches above or below the graph, if 
this is appropriate. This option can be used only if 
the Y-scaling is automatic (YSCALE = 0). 
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= 2 


The X-axis location will be as specified by IXUP. 


m. MODYAX: 


Determines the mode of Y-axis location in the same 
way as MODXAX, above, governs the X-axis location. 


n. IWIDE: 


Width of graph in inches (1^IWIDE<9). If IWIDE is 
out of this range, a value of 8 will be assumed. 


o. IHIGH : 


Height of graph in inches (IS IHIGH S 15) . If IHIGH 
is out of range, a value of 8 will be assumed. 


p . I GRID : 


If IGRID = 1, a 1" x l 11 grid will be superimposed on 
the graph. This is useful only if plain paper is 
used on the CalComp Plotter. 


q . LAST : 


Indicates to the calling program whether the previous 
plot was completed successfully. The codes are: 


= 0 


Last plot was completed successfully 


= 1 


Last plot was not completed successfully. 


= 2 


Last plot was not completed successfully and no further 
graph output will be attempted until DRAW is next 
entered with MODCUR = 1 or 0 


= 3 


An attempt was made to enter DRAW With MODCUR f 1 or 0 
while the error lockout was set. This argument must 
always be a variable name and never a number in the 
call statement. 
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3. Notes and Ccmments : 



a. The graph scales and, if MODXAX = 1 and/or MODYAX = 1, the 
amounts of origin offset are always output as part of the 
graph title. 

b. There are internal checks of the input to DRAW to prevent 
incorrect use. If an input error is detected, an attempt 
will be made, where possible, to complete the plot. If 
an argument is "correc ted M in this process, the user will 

be so informed on the standard output. If it is not possible 
to complete the plot, the user will be informed of the 
reason by a message on the standard output. 

c. If part or all of a curve would fall more than 0.6" laterally 
beyond the ends of the X-axis , or 0.7" vertically beyond 

the ends of the Y-axis, the X and/or Y ordinates will be 
limited so that the curve will typically become a line along 
part or all of the boundary of the graph as here defined. 

d. If one or more points of a point plot would fall outside the 
graph area, the plot of the point, or points will be inhibited. 
The number of such points will be reported to the user on the 
standard output. 

e. It should be pointed out that the X and Y scaling and axis 
locating processes are entirely independent, so that, for 
example, X might be auto-scaled, while the Y-scale is 
specified. At the same time the X-axis might be located 
automatically, while the Y-axis location is specified. 

f. It must be remembered that the scales and axis locations of 
a multi-plot graph are set when DRAW is called for the first 
time (with MODCUR. = 1). Thus the user must attempt, at the 
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time to achieve scaling and axis location which will be 
appropriate to all the plots he intends to make on the 
one graph. Particularly if the automatic features of 
DRAW are selected, foresight will be demanded of the user 
in this respect. 

4. Auto- Scale Properties: 

The scale factor is chosen from amongst the values 1, 2, or 5 
units per inch, or some power of 10 times one of those values. 

A curve, or set of points which is plotted with auto-scale will 
normally lie entirely within the graph area is defined in 3.d., 
above. The only exception may occur if an axis is placed, by 
the user, along one edge of the graph (e.g., IXUP = 0, MODXAX = 2). 
In such a case, points "outside 11 the axis are not considered in 
the selection of a scale factor (e.g., negative do not affect 
the choice of scale when IXUP =0). If automatic axis location 
as well as auto-scale is selected, the plot, if it does not fill 
the graph area, will be placed as far as possible towards the 
bottom- left of the graph area consistent with the fact that the 
axes can be set in increments of 1" only. 

1st revision - 31 July 1967. 
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APPENDIX B 

COMPUTER PROGRAM FLOW CHART 
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END 
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APPENDIX C 



DATA CARD DESCRIPTION 
C-l Data Card Placement 

For one set of data, six data cards are required to give the 
program the necessary input. The first card gives the physical data 
such as leg length, pipe size and material, and other quantities that 
describe the physical situation. The second card supplies the control 
variables to the program. These control variables indicate to the 
program what parameter is to be investigated and its range. 

The last four cards supply the initial values given to the 
FTE f s (see page 26). Each of the last four cards contain eight of 
these values. Unless a particular computer run starts near the 
critical part of the parameter range (where convergence is difficult) 
the values on the last four cards are usually zero. If such is the 
case, blank cards may be used; however, when values other than zero 
are supplied as input, the format shown in Appendix C-2 must be used. 

It must be remembered that six data cards must always be used 
for a particular set of input whether or not the FTE ? s are zero. 

Any number of sets of six data cards can be supplied to the 
program for any one run. Consideration should be given to time 
requirements, however. After the final set of six data cards, one 
blank card follows and then, after that, the last control card (see 
Appendix D) . The blank card is used to terminate the run by telling 
the program that all data has been considered. 

The user should be familiar with fixed and floating point variables 

3 

and the common Fortran field specifications. 

3 Daniel D. McCracken, A Guide to Fortran IV Programming (New York: 

John Wiley and Sons, Inc., 1965), pp. 3-5, 13, 84. 
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Appendix C-2 shows a typical input (all six data cards). The 
extreme left hand column indicates the data card number (this order 
must always be followed). For each data card there are two rows 
shown. The top row indicates the variable name whose value is 
shown directly below in the next row. This second row is the 
actual value of the variable named above it and appears just as the 
data card itself does. The top row of the two in each case serves 
only to aid in identification of the variable. 
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C-2 



Sample Data Card Set 
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Blank 

one SET 

OF 

SIX 

DATA 

CAR.OS 




APPENDIX D 

CONTROL CARD ARRANGEMENT 
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APPENDIX E 
PROGRAM LISTING 



E-l Main Program 
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E-2 Subroutine ALPHA 
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E-3 Subroutine FACTOR 
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CIN OOUT DTEMP P E COEFI K PAR 1 
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ELAPSED TIME = 16. C5 



APPENDIX F 



SAMPLE CALCULATION 



The following table of values is provided to aid the user in 
checking the computer calculations. Two combinations of secondary 
effects are listed. The first is the "Exact Value" combination 
where all secondary effects are included. The second is the most 
erroneous combination, or the one where all secondary effects are 
suppressed in the calculation. The descriptive data is that listed 
on page 42 of Section IV and also exactly the same as the sample 



data card in Appendix 
Quantity 
Assumed FI 
Assumed F2 
W 

<*S 

<Xio 

WBAR 

5c, 



C, with a leg ratio of 
All Effects Allowed 
65372.66553 
4623.212071 
.0054419338 
.5011360506 
.3340605064 
1.001817935 
1.002727150 
.4185170188 
.2678453147 
.0072796731 
.6695786688 
.4320902252 
.2859409726 
1.923736751 
2.324235399 
1.496163752 



. 2 . 

All Effects Suppressed 
74373.10291 
5192.148178 
0.0 

.5000000000 
.3333333135 
1.000000000 
1.000000000 
.4166666865 
.2666666508 
0.0 

.6666666865 
.4166666865 
.2666666508 
0.0 

.5000000000 
.3333333135 



89 



Quantity 



All Effects Allowed 



*3 

5*4 

&S 

5ti 

o( e 

0< (0 

*1 

< 3 ,, 

Q\Z 

^21 

^at 

A 

A 

FI Converged 
F2 Converged 
Moment 



3.878225194 

5.471216862 

8.613517662 

5.485230450 

31.81243915 

13.54446691 

8.624985666 

5.499880474 

1.817163214 

205.4115482 

-.0001242864 

.0018232014 

.0000294904 

-.0001211686 

.0132454070 

.0132186939 

65840.74973 

4633.905400 

1071501.071 



All Effects Suppressed 
1.000000000 
1.000000000 
.4166666865 
.2666666508 
0.0 

.6666666865 

.4166666865 

.2666666508 

.6000000000 

150.0000000 

-.0000882817 

.0013242253 

.0000247189 

-.8828170315 

0.0 

0.0 

74972.17486 

5208.593999 

1231122.149 
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APPENDIX G 



SHORTENING DUE TO DIFFERENCE BETWEEN ARC AND CHORD 
The secondary effect produced because of difference between 
arc and chord is described by equation (7) of Section II when 
higher order terms are neglected. The following is a more detailed 
development of that relationship. 

In order to have a physical understanding of the background of 
the derivation, refer to Fig. G-l. Due to the action of cross shear, 
as it is idealized in customary theory, the horizontal length, dx, 
from 0 to A, is stretched as it is inclined to angle , so that it goes 
from 0 to B. Then the action of bending further inclines it so that 
it takes the position OC, but there is no change in length during this 
bending process. (Change in length due to axial force is accounted 
for separately.) The increment of difference between arc and chord 
is the distance dA shown in the figure. 

In the derivation, it is first assumed that Y is small so that 
tantf approximately. 




d s z = dx 2 -H dtj 1 

dA * dx - ds cos(^) 

^ = * or dy - tfdX 

ds* = dx* +■ 

ds « dx* ( i + X*) 

♦\ ds » dx fi"+ 7 r 

d>\- dx-dxfi+V* Cos(y') 

dh- dx()-fi+t T CosCs!)) 
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is substituted 



If series representation of cos(y ! ) and \ IT+\ r 4 
with higher order terms neglected 

dX= 1 IT (('if -t) 

or 

^ = i ^ [(a 1 )* - - * l ]dx 



( 7 ) 
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